Myosin V passing over Arp2/3 junctions: branching ratio calculated from the elastic 

lever arm model 
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Myosin V is a two-headed processive motor protein that walks in a hand-over-hand fashion along 
actin filaments. When it encounters a filament branch, formed by the Arp2/3 complex, it can either 
stay on the straight mother filament, or switch to the daughter filament. We study both probabilities 
using the elastic lever arm model for myosin V. We calculate the shapes and bending energies of 
all relevant configurations in which the trail head is bound to the actin filament before Arp2/3 and 
the lead head is bound either to the mother or to the daughter filament. Based on the assumption 
that the probability for a head to bind to a certain actin subunit is proportional to the Boltzmann 
factor obtained from the elastic energy, we calculate the mother/daughter filament branching ratio. 
Our model predicts a value of 27% for the daughter and 73% for the mother filament. This result 
is in good agreement with recent experimental data. 



I. INTRODUCTION 

Myosin V is a two-headed processive motor protein 
from the myosin superfamily, involved in different forms 
of intracellular transport [H, 0- It has drawn a lot of 
attention in recent years and is now one of the best stud- 
ied motor proteins. The experiments have characterized 
it mechanically 

i, II 0, MM E, biochemically 
11. [l2 L Ha . fl4i . optically [15|, [l(| LL21 and structurally 
la . [l9l I2Q , l2l| . These studies have shown that myosin 
V walks along actin filaments in a hand- over- hand fash- 
ion [TBI, Hi] with an average step size of about 35 nm, 
roughly correspon ding to the periodicity of actin fila- 
ments d i, H, [H, Ei, a stall force of around 2 pN Q 
and a run length of a few microns [1, Hi, [24|. Under 
physiological conditions, ADP release has been identified 
as the time limiting step in the duty cycle [H, HH • 

The Arp2/3 complex HH] initiates actin filament 
branching in the vicinity of a protruding edge of a 
cell. The complex consists of 7 subunits (Arp2, Arp3 
and ARPC1 through ARPC5) and is activated by 
WASp/Scar proteins Hi]- It binds to the side of one 
("mother") actin filament and initiates the nucleation 
of a second ("daughter" filament), which starts grow- 
ing with the fast growing end ("+" end) away from the 
Arp2/3 complex. The mother and the daughter filament 
enclose an angle of 70° . 

The question what happens to a myosin V motor when 
it arrives at an Arp2/3 mediated actin filament junction 
is of interest for several reasons. First, the branching 
behavior is important for understanding vesicle transport 
in the actin cortex. And second, it is of high interest 
when studying the fundamental mechanism of myosin V 
stepping, because it represents a well defined situation in 
which predictions from different theoretical models can 
be tested against the experimental result. Our aim in this 
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article is to use the elastic lever arm model for myosin 
V, which is described in detail in a previous article [27l ]. 
to predict the dynamics of a myosin V motor that passes 
over an Arp2/3 junction. In particular, we will calculate 
the probabilities that a motor continues along the mother 
or the daughter filament. 



II. MODEL 

The idea behind the elastic lever arm model for myosin 
V is to describe the dimeric motor as an assembly of two 
identical heads, connected together and to the cargo- 
binding tail with elastic lever arms. The model allows 
us to derive the properties of a dimeric molecule, such 
as step size distribution, force-velocity relation and pro- 
cessivity, from the properties of an individual head, such 
as geometry, chemical kinetics and elasticity. In this re- 
spect the approach is different from the class of discrete 
stochastic models, which describe the motor as a single 
unit HI HI]. 

We describe each head with a 5-state mechano- 
chemical model, similar to that for muscle myosin (e.g., 
Hi]), where each state (with bound ADP.Pi, ADP (pre- 
powerstroke) , ADP (post-powerstroke), without a nu- 
cleotide, detached) is connected with a certain orienta- 
tion of the lever arm, as determined with electron mi- 
croscopy (EM) [lEGj]. Th e lever arms are modeled as 
elastic beams, connected with a flexible joint. A recent 
study measuring fluctuations in the position of the free 
head H3 has demonstrated that the lead head diffuses 
around the joint freely before binding to the next actin 
site, meaning that there is no detectable elastic energy 
cost connected with variation in the angle between the 
two lever arms. The very nature of protein flexibility, 
which mainly originates from the twisting of bonds be- 
tween carbon atoms in the backbone, leads us to the 
conclusion that the joint is also fully flexible with regard 
to rotation of each lever arm along its axis. Similar flex- 
ibility has also been observed in myosin II H3, Hi, H3 • 

The calculation of the branching probability is simpli- 
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TABLE I: Geometric parameters of the Arp2/3 junction and 
a myosin V head. 



Distance actin subunits 


a 




2.75 nm 


Angle actin subunits 


o 




167.14° 


Daughter filament angle 


P 




70° 


Daughter filament offset 


r 




12 nm 


Daughter filament rotation 


7 




39° 


Lever arm start: radial pos. 


R 




8 nm 


Lever arm start: displacement 


^pre — 


PS 







5post 


-PS 


3.5 nm 


Lever arm angle 


0pre- 


PS 


115° 




<Ppost 


-PS 


50° 


Lever arm length 


L 




26 nm 



tied a lot if we make the following assumptions. First, we 
assume that the binding of the lead head always leads 
to a step, which means that its unbinding is significantly 
slower than the step that follows in the regular cycle (Pi 
release). Second, we assume that the probability for the 
lead head to bind to site j if the trail head is bound to 
site i is given by the Boltzmann factor 



Pa 



exp 



k B T 



(1) 



Ui.j denotes the elastic energy of deformed lever arms 
when the trail head is in the post-powerstroke state, 
bound to site i, and the lead head in the pre-powerstroke 
state to site j. In the following, we will use the notation 
where the sites on the mother filament are marked with 
(M, i) and those on the daughter filament with (D,i). 
For example, Pd,2|m,-9 denotes the conditional proba- 
bility for the lead head to bind to site 2 on the daughter 
filament if the trail head is bound to the site —9 on the 
mother filament. We enumerate the actin subunits so 
that the central subunit under the Arp2/3 complex on 
the mother filament has the index 0. Positive indices 
denote subunits towards the "+" end and negative to- 
wards the "-" end. Subunits of the daughter filament are 
enumerated from onwards. Note that sites numbered 
-2, and 2 on the mother filament are not accessible for 
a myosin V head, because of steric hindrance with the 
Arp2/3 complex. 

The structure of the Arp2/3 complex and both actin 
filaments (Fig. [I} has been determined from EM stud- 
ies [35l |3(1 ] and its parameters are summarized in Table 
HJ While we approximated actin with the commonly as- 
sumed 13/6 helix in the original paper [27|, we use a more 
accurate 28/13 helix here, with the angle do = 167.14° 
between adjacent subunits. A detailed discussion on dif- 
ferent helix models and their consequence for the calcu- 
lated step size distribution can be found in [3?| . 

When a head is bound to the site i, the starting point 
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FIG. 1: Arp2/3 junction and a dimeric myosin V motor. (A) 
The Arp2/3 complex is attached to the side of the mother fil- 
ament occupying subunits -2, and 2. It nucleates the growth 
of a daughter filament, whose position is determined by the 
angles (3 (branching angle) , 7 (twist) and r (distance from the 
first subunit to the center of the mother filament). (B) The 
myosin V motor is consists of two heads, connected with lever 
arms, which we describe as elastic rods. The proximal end of 
each lever arm always leaves the head at a fixed angle <f>, which 
depends on the nucleotide state of that head. The distal ends 
of both lever arms are connected with a fully flexible joint. 



and the unit vector giving the initial direction of the lever 
arm are given by 



x = R x (-i8 ) 



'ia + A /cos< 

t o = B*(-»0o) I (2) 
\ R J \sin 



where R x denotes the rotation matrix around the x axis, 



R*(0) 



/ 1 

cos 8 — sin ( 
V sin 9 cos 8 



(3) 



For a head bound to the daughter filament, the two 
vectors read 

Ir + ia + 5\ 
xo = R„(-0)Rx(7-*0o) 



to = Rj,(-/3)R x (7-^o) 



cos 




(4) 
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with 

(cos (3 sin (3 \ 

10. (5) 

— sin cos (3 J 

Here (3 = 70° denotes the angle between the mother and 
the daughter filament and 7 = 39° the rotation of the 
daughter filament around its axis (see Fig. [IJ . 

We calculate the shapes of both lever arms as de- 
scribed in Ref. (27j by minimizing the bending energy 
U = J ds EI(C(s)) 2 /2, where C(s) denotes local curva- 
ture. For the bending modulus of the lever arm we use 
the value EI = 1500 pNnm 2 , which corresponds to a 
"spring constant" of k = 3EI/L 3 = 0.25pN/nm, mea- 
sured at the tip of the lever arm. This value was orig- 
inally estimated from the stall force, but it shows good 
agreement with direct optical tweezer measurements [38| . 
We neglect any additional compliance resulting from the 
head or converter domain. Because most of the bending 
takes place in the proximal part of the lever arm, we ex- 
pect that its effect would not be significantly different. 
For the temperature we use the value T = 27°C. 

III. RESULTS 



TABLE II: Probability distribution Pi,i for the first accessed 
binding site in the interval —17 < i < —3. Pj,i denotes 
the probability distribution for the first accessed site in the 
interval —13 < i < 1. Pj^ is calculated according to Eq. (|ll|l . 
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0.0759 
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(Pl3 + Pl5 
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0.0758 


-2 








-1 




(Pis + Pn)/l 


0.0754 











1 




Pis/r 


0.0075 



Stepping of myosin V in the absence of the Arp2/3 
complex 

When a myosin V motor is sufficiently far away from 
the Arp2/3 complex, it exhibits the step ping pattern that 
has already been discussed in Refs. [23, [33]. We restrict 
the step lengths of an unperturbed motor to 11, 13 and 
15 subunits, with probabilities 

P11 = Pu,i+ii\M,i » 0.005 , (6) 
P13 = Pu,i+i3\M,i ~ 0.895 , (7) 

-Pl5 = fM,i+15|M,i ~ 0.1 , (8) 

calculated from Eq. |T]). Probabilities for other step sizes, 
such as Pg and Pyj turn out to be very small. Note that 
in this calculation P\\ and P15 are likely to be somewhat 
underestimated - their values are somewhat higher if we 
take into account torsional fluctuations in the actin helix 

The average step size can be calculated from these 
probabilities as 

f=llPii + 13Pi3 + 15Pi 5 «13.2. (9) 

In the absence of the Arp2/3 junction, the fraction of 
sites that get accessed by a passing myosin V motor is 1 /I. 
This is also the probability that a site before the junction 
(i < —7, as will be shown later) ever gets accessed by the 
motor: 

P M a = j for % < -7 . (10) 



Initial state: First accessed site in the interval 

-17 < i < -3 

We start our analysis at the point where a head of an 
approaching myosin V first passes the subunit —17 or 
binds to it. This way the initial state is definitely not 
influenced by the presence of the Arp2/3 complex. By 
counting only the first accessed site in this interval, we 
avoid double-counting of events where the motor binds, 
for example, first to site -17 and then to -4. 

We denote with Pjj the probability that once the first 
head has bound to any subunit in the interval — 17 < i < 
—3, this has happened at subunit i. P\.i can be calculated 
in a way that is illustrated in Fig. [5] The sites between 
-17 and -7 can only be reached from outside the interval. 
Therefore, whenever the motor binds to one of them, 
it becomes the first accessed site in this interval. The 
probability Pij then equals the probability that the site 
i ever gets accessed by the motor, which is l/l, according 
to Eq. pTJ)) . For i = —6 and higher the situation becomes 
different. Site -6 counts as the first site in this interval 
if the preceding step size is 13 or 15 subunits, but not 
if it is 11. Therefore, the corresponding probability is 
Pi,_6 = (P13 + Pin)/l- The probability that the first site 
accessed in this interval is -5, Pi,-5, has the same value. 
Finally, the sites -4 and -3 will be the first accessed sites 
in the interval if they are following a 15 subunit step. 
Therefore, their probabilities are Pi,_4 = Pi. -3 = P15/L 
The probabilities P\^ are given in the second column of 
Table UT1 and shown in the top graph of Fig. [31 
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FIG. 2: Different pathways on which a myosin V motor can approach the Arp2/3 junction. The probability distributions for 
the first accessed site in the two intervals marked with dashed lines are denoted as Pij and Pj t i. Pi t i represents the probability 
that i is the first accessed site with — 17<i<— 3 (between dashed lines in the first or second column). For values between 
—17 and —7, this state can be reached with 3 different step sizes. For i — —6 and i — — 5, it can be reached with 2 different 
step sizes (13 and 15). If it is reached by a shorter step (11 subunits), it means that the preceding binding site was already 
inside the interval —17 ... — 3, so i is not the first accessed site within it. For i — — 4 and i = — 3, the site can only be reached 
with a 15-subunit step. Pj^ denotes the probability that i is the first accessed subunit in the interval — 13 < i < 1 (dashed 
lines in the right column). 




FIG. 3: Probability distribution Px : i (upper diagram) for the 
first accessed binding site on the mother filament in the inter- 
val — 17 < i < —3. The lower diagram shows the probability 
distribution Pj,, for the first accessed site with —13 < i < 1. 
Pj ; j is calculated according to Eq. (|llfl . by redistributing the 
probabilities for —17 to —14 to other sites, as indicated by 
arrows. For example, if -17 is the first accessed site with 
—17 < i < —3, the first accessed site with — 13 < i < 1 can 
either be -6 or -4. Note that the total probability that the 
motor binds to site -4 is higher than for any other site, which 
is due to the inaccessibility of site -2. 
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FIG. 4: Most probable lead head binding sites for different trail head positions. For selected trail head binding sites i, a 
collection of lead head binding sites and their bending energies Uij (in units of fcsT) is shown. Some trail head positions which 
do not have significant branching probabilities (e.g., —12) are omitted. The shape of both lever arms in each configuration is 
calculated numerically, by minimizing the bending energy. 



As a next step, we will derive the probability dis- 
tribution Pj t i for the first accessed site in the interval 
— 13 < i < 1. We choose this interval because (M,-13) is 
the first site from which the motor can reach the daugh- 
ter filament. This distribution can be obtained from P\ t i 
by redistributing probabilities for sites between -17 and 
-14 according to the conditional probability for the next 
step 



P 



3,3 



*j + E 



PM,j\M,iPl,: 



(11) 



i=-17 



as shown in Fig. [3) The values of Pjj are given in the 
third and fourth column of Table|TTJ This distribution de- 
fines the state from which we will calculate the branching 



ratio at the Arp2/3 junction in the following section. 



Conditional branching ratio 

Now we can calculate the conditional probabilities that 
the lead head binds to the daughter filament, if the trail 
head is bound to a mother filament subunit i fulfilling 
— 13 < i < 1. We denote this probability as Pd\M,i- 

Figure 0] shows the most relevant dimer configurations 
with the trail head on mother filament sites between -15 
and 6 and the lead head either on mother or on daughter 
filament. For each trail head position, the conditional 
probability that the lead head binds to a certain site is 
given by Eq. ([T]) with the index j' running over all ac- 
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FIG. 5: An example of a one-step (upper path) and two-step (lower path) process through which the myosin V motor can 
reach the daughter filament. The probability that the trail head is bound to subunit (M,-9) in the initial state is Pj,-g. In the 
next step, the lead head can bind to site (D,2) with conditional probability -Pd,2|m,-9- We denote such processes as one-step. 
Alternatively, the lead head can bind to the site (M,4) with conditional probability Pm,4|m,-9 and, subsequently, the other head 
can bind to site (D,8) with conditional probability Pd,8|m,4- This is an example of a two-step process. Note that this scheme 
shows just two examples of possible pathways and omits alternatives that are indicated by dashed arrows. 



cessible mother- as well as daughter-filament sites - cor- 
responding to one row in Fig. [4] Generally speaking, 
the daughter filament can either be reached directly, in a 
one-step process such as (M,-6)^(D,3), or in a two-step 
process such as (M,-9)^(M,4)^(D,8). Two examples 
are shown in Fig. O The probability that a motor with 
the trail head bound to site (M,i) eventually binds to the 
daughter filament is 



Pd\M,i — 



j 



P 



DjjM, 



M,k\M,i 



(12) 



Here the first term denotes one-step processes and the 
second term two-step processes, where the motor first 
binds to site k on the mother filament and subsequently 
to site j on daughter filament. Our numerical calcula- 
tion shows that the only significant terms are those with 
k — 4. The branching ratio Pd\yi,i f° r each trail head 
position is shown in Fig. [5] The graph shows separately 
the contributions of one- and two-step processes. 
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FIG. 6: Probability that the motor will step on and continue 
its walk along the daughter filament if the trail head is ini- 
tially bound to site i on the mother filament, Pd\yi,i (black 
and hatched bars), as calculated from Eq. (|12p. The black 
bars show the contribution of one-step processes, in which 
the lead head binds to a site on the daughter filament imme- 
diately. The hatched bars show the contribution of two-step 
processes in which the motor first binds to another site on 
the mother filament (usually 4) and then in the second step 
to the daughter filament. The grey bars show the probability 
that the motor continues along the mother filament. 



Total branching ratio 



IV. DISCUSSION 



With these probabilities and weights Pj^ we finally 
obtain the branching ratio for the daughter filament: 



Ph = 



i=-13 



Pj,iPd\M,i 



0.27 



(13) 



If one head binds to the daughter filament, there is still 
some probability that the next head will bind back to the 
mother filament. One such example, with the trail head 
on the site (D,l), is shown in the last row in Fig.[U How- 
ever, the contribution of such events to the total branch- 
ing ratio is not significant. 



While the exact result calculated above does require 
to take into account all the individual configurations, 
its order of magnitude can also be understood with a 
simple "handwaving" argument, which goes as follows. 
Roughly speaking, the approaching myosin V motor can 
either reach the Arp2/3 complex on the opposite side of 
the actin filament, in which case it cannot switch to the 
daughter filament, or on the same side, in which case the 
probability to switch to the daughter filament is about 
1/2. Together, this gives a branching ratio of 1/4, not 
far from the exact result of 0.27. 

An experiment measuring (among other quantities) the 
branching ratio at Arp2/3 mediated actin filament junc- 
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tions was recently carried out by Ali et al. [39|. The 
results are not directly comparable - in the experiment 
the actin filaments were attached to a glass surface so 
that some binding sites were not accessible for myosin V. 
However, because the blocked sites are different depend- 
ing whether the side filament branches to the left or to the 
right, we expect that in statistical average, the calculated 
branching ratio still gives a good approximation. In the 
experiment 18% of the molecules dissociated, 20% con- 
tinued on the daughter filament, and 62% on the mother 
filament. If we discount dissociation events, this means 
that the fraction that switched to the daughter filament 
was 24%. The statistical error of this figure is about ±5% 
(the total number of events observed was 76). Therefore, 
the result can be considered in excellent agreement with 
the model calculation. 

Note that our calculation only concerns stepping pat- 
terns and does not take into account kinetics. This has 
an important advantage that the result only depends on 
geometric and elastic parameters, but not on kinetic con- 
stants, some of which are still less well known [27j . Com- 
bining the present model with the full kinetic scheme 
could influence the predicted branching ratio as follows. 
First, the binding rate of the lead head can be different 
in the vicinity of the Arp2/3 complex. And second, the 
ADP release rate of the rear head can be influenced by in- 
tramolecular strain, possibly also by its lateral (off-axis) 
component, as observed by Purcell and coworkers [4(ij |. 
Both these effects do not have a direct influence on the 
branching ratio, but they might have an indirect one by 
influencing the dissociation probability. In the present 
calculation, events where the whole myosin V molecule 
dissociates from the actin filament are not taken into ac- 
count. There is, however, evidence that the predominant 
dissociation path leads through detachment of a head in 
the ADP state j23l . 1 4 ill - these processes are denoted as 
Pathway 2 in Ref. |27| . It is therefore plausible that the 
termination rate increases if either the lead head is hin- 
dered in its search for a binding site, or the ADP release 
in the trail head is slowed down. In our model (see elastic 
energies in Fig. [4J , the lead head binding rate is strongly 
reduced if the trail head is bound to site —13, it is also re- 
duced somewhat if the trail head is bound to —11, while 
it is accelerated for trail head positions —9 and 4. In 
total, dissociation can be accelerated by the presence of 
the Arp2/3 junction in 2 out of I cases. Without going 
into quantitative details, one can conclude that the dis- 
sociation probability could theoretically increase by up 
to 2/1 « 0.15. The branching ratio for the daughter fil- 
ament could then be somewhat smaller, because those 
trail head positions that have the higher branching ratio 
are also more likely to lead to dissociation. The effect 
caused by the strain-dependent ADP release rate is more 
difficult to estimate, mainly because the exact influence 
of off-axis strain on the ADP release is not yet known 
quantitatively 

To check the robustness of our result against uncertain- 
ties in model parameters, we calculated the variation of 
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FIG. 7: Dependence of the calculated branching ratio Pd on 
different model parameters. (A) Dependence on the lead head 
angle <^ pre -ps (solid line, lower scale) and trail head angle 
0post-ps (dashed line, upper scale). (B) Dependence on the 
lever arm stiffness (bending modulus) EI. In both diagrams, 
the dot shows the value used in all other calculations through- 
out the article. 



the branching ratio with several model parameters. Note 
that these calculations were carried out numerically and 
took into account all relevant processes, including longer 
and shorter steps, as well as transitions not shown in 
Fig. 01 Data in Fig.[7]shows a variation of about ±10% if 
either the lead head or the trail head angle is modified by 
±15° . The allowed variation of these two parameters that 
keeps the model consistent with the experimental result is 
therefore restricted, although parameter sets where both 
angles are increased or decreased simultaneously cannot 
be excluded. The result is more robust against variations 
in the lever arm stiffness EI, where deviations do not ex- 
ceed few percent if the value of EI is changed by a factor 
of 3 in either direction. 

We can therefore conclude that the calculated branch- 
ing ratio adds support to the elastic lever arm model 
presented in Ref. [27] and the geometric parameters used 
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there. However, we cannot use it as a criterion to de- 
termine the lever arm stiffness, which is still not known 
precisely. Another open question is to what extent the 
result can be reproduced with alternative models, such as 
Ref. [42|], which uses a more complex model for the elas- 
ticity of the lever arm, with a soft longitudinal (about 1/3 
of the value used there), but very stiff azimuthal compo- 
nent. The completely different class of "hot spot" mod- 
els, which proposes that the position of the next bind- 
ing site is determined by a propagating conformational 
change in the actin filament [431 ] . on the other hand, 
seems less compatible with the finding, unless the con- 



formational change could propagate through the Arp2 / 3 
complex as well. 
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